Quantum cavitation in liquid ^He: dissipation effects 



o\ 



B 

a 
o 
o 



1 



Dora M. Jezek^, Marti Pi^, and Manuel Barranco^. 
^Departamento de Fisica, Facultad de Ciencias Exactas y Naturales, 

Universidad de Buenos Aires, RA-1428 Buenos Aires, and 

Consejo Nacional de Investigaciones Cientificas y Tecnicas, Argentine 

"^Departament d'Estructura i Constituents de la Materia, Facultat de Fisica, 

Universitat de Barcelona, E- 08 028 Barcelona, Spain 

Q\ • (February 1, 2008) 

0^■ 



§ ■ Abstract 



We have investigated the effect that dissipation may have on the cavita- 

J^ I tion process in normal liquid ^He. Our results indicate that a rather small 

C^ ' dissipation decreases sizeably the quantum-to-thermal crossover temperature 

^N ■ T* for cavitation in normal liquid ^He. This is a possible explanation why 

f~^ . recent experiments have not yet found clear evidence of quantum cavitation 

On , at temperatures below the T* predicted by calculations which neglect dissi- 
0\ 



pation. 
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Quantum cavitation in superfluid liquid ^He has been unambiguously observed using 
ultrasound experimental techniquesilB. These experiments have shown that quantum cavi- 
tation takes over thermal cavitation at a temperature (T) around 200 mK, in good agreement 
with theoretical calculationsEfcl, so that the problem of cavitation in liquid ^He can be con- 
sidered as satisfactorily settled. 

The crossover temperature corresponding to ^He has also been calculated&u, predicting 
that T* ~ 120 mK. It turns out that preliminary results obtained in a recent experimental 
have not shown clear evidence of quantum cavitation for temperatures even below that value. 
However, the phenomenon has been firmly established as a stochastic process. A possible 
explanation is that thermal cavitation is still the dominant process down to temperatures 
lower than predicted. 

The method of Ref. |^ (see also Ref. ||) is based, on the one hand, in using a density 
functional that reproduces the thermodynamical properties of liquid ^He at zero temperature 
(equation of state, effective mass, etc), as well as the properties of the ^He free surface. A 
major advantage of using a density functional is that one can handle bubbles in the vicinity 
of the spinodal region, where they are not empty objectsSa and any attempt to describe the 
critical bubble in terms of a sharp surface radius failsQ. On the other hand, we have used a 
functional- integral approach especially well suited to find T*. This gives us some confidence 
on the values obtained for the crossover temperature, and inclines us to think that any 
appreciable discrepancy between theory and experiment has to be attributed not to the 
method itself, but to some physical ingredient which has been overlooked in the formalism. 
One such ingredient in the case of liquid ^He is dissipation, which is knownu to decrease 
T*. Since ^He is superfluid below the lambda temperature, we are actually treating both 
quantum fluids within the same framework, the behavior of ^He being accounted for by the 
dissipationless version of the general formalism. 

Our starting point is the real time Lagrangian density C{p, s) 

C{p,s) = mps-n{p,s) , (1) 

where p{r,t) denotes the particle density, m the ^He atomic mass, and s{r,t) is the velocity 
potential, i.e, the collective velocity is u{r,t) = Vs(r, t). The Hamiltonian density 7i(p, s) 
reads 

n{p,s) = -mpu^ + [uj{p)-uj{pm)] , (2) 

where uj{p) is the grand potential density of the system and pm is the density of the 
metastable homogeneous liquid. We refer the reader to Ref. ^ and references therein for 
details. 

To describe the dynamics in the dissipative regime while still being able to deal with 
inhomogeneous ^He, which is crucial for a proper description of cavitation in liquid helium, 
we have introduced a phenomenologycal Rayleigh's dissipation functionpOJ JF 
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From Lagrange's equations 
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with X being either s or p, one gets the continuity and motion equation, respectively: 

p + V{pu) = 



(4) 



(5) 



jduk 
m < -^ + UiVkUi 



-^^0+^^ 



^;^v(p«) 
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For an homegeneous fluid, the equation of motion ressembles the Navier-Stokes equation! 







f duk ^ 
mp< — — + UiVkUi 



-VkP + V^Uk + ( C + 3 ^ ) Vfc(V ■ u) 



(7) 



where P is the pressure. For hquid ^He at low T, dissipation depends on the mean free path 
of quasiparticles, and a precise estimation of the magnitude of this effect in the tunneling 
process is difficult. Since our interest here is to explore the effect of a small viscosity on T*, we 
have adopted the pragmatic point of view of identifying ^ with ^ + 77/8 and presenting results 
for different ^'s close to the experimental r] value (it is known that at low temperatures, the 
shear viscosity coefficient r] is much larger than the bulk viscosity coefficient C,., see for 



example Ref. 12). Using the macroscopic viscosity coefficient, one should have in mind that 
we are likely overestimating the dissipation effects. 

To obtain T* we have proceeded as indicated in Ref. ^, writing the above equations in 
imaginary time t = it and linearizing them around the critical bubble density pq, seeking 
solutions of the kind: 

p(f,r)=Po(r) + p^(r)e-^-»^. (8) 

Upon linearization, we end up with the following equation for Ug and p^{r): 

Mp\r) = [muJ^^-Mi~ ^uJsM2\ p\r) = . (9) 

The differential operators 7VI1 and A^2 in Eq. ® are, respectively, the linearization of 

and V l^pV (^1) I , (10) 

in which only first order terms in p^ (r) and its derivatives have been kepw Since ^ depends 
on the densityliHl as p^^^, in actual calculations we have made a local density approximation, 
using as form factor in Eq. (^) the expression I /{PsLp^^^ {''")) ^ where psat is the density of 
the liquid at T = and P = 0, and ^ is then density-independent. 

Eq. (^ is a fourth-order linear differential, generalized eigenvalue equation, whose phys- 
ical solutions have to fulflU p^'(O) = p^'"(0) = 0, and fall exponentially to zero at large dis- 
tances. We have solved it as indicated in Ref. ^. Once the largest dissipation-renormalized 
frequency tUs has been determinedlid, the crossover temperature is obtained as T* = ?i(jJs/{2tc). 




Table I collects the equation of state near the spinodal point {psp = 0.01191 A , Psp = 
—3.102 bar), and other quantities which are of interest to analyze the experimental resultso. 
Our spinodal point compares very well with recent Monte Carlo calculationalil {psp = 0.0121 
A , Pgp = —3.12 ± 0.10 bar), and also with other phenomenological approached^EJ. 

We show T* in Fig. |l| as a function of pressure for different ^ values. In particular, ^ = 
100 fiP roughly corresponds to the experimental valuecB of ?7 at P = and T = 100 niK. The 
associated effective quantum action S obtained as 5 = AQ/T*, where AQ is the maximum 
of the energy barrier, is displayed in Fig. 0. 

Fig. ^ shows p^(r) at P = — 3 bar for three ^ values, as well as the critical bubble 
density po{r). The linearized continuity equation p^(r) oc V(po^) implies that p^(r) must 
have nodes, as it imposes that the integral of p^ (r) is zero when taken over the whole space. 

When ^ is small enough and the A^2 term in Eq. (j^) can be treated perturbatively, a 
straightforward calculation yields 



Col., 



2 , f^l^A ^/i2 .....v 



where we have used a standard matrix notation&Ej to denote as u;o,o and \pq ) the higher 
frequency solution of the non-viscous problem {raujQ^ — M.i)\p\^^^) = 0, and have defined 

P2 = — (Po 1-^2 IPo ) > 0. Equation (|TT]) is similar to that given in Ref. | for the 
dissipation-renormalized frequency ujg in the case of frequency-independent damping. 

Figures (|I|-0) indicate that for viscosity values of the order of the experimental one, 
a sizeable decrease of the crossover temperature occurs. However, the present model still 
predicts that a transition from thermal to quantum cavitation takes place in liquid ^He. 

We finally obtain the homogeneous cavitation pressure Ph from the equationQ: 

1 = (W)e.p Jo e-'' , (12) 

taking for the experimental volumextime {Vt)exp a typical value of 10^ A s, which corre- 
spods to "^He experimentsD. We have adopted for Jq the same prescription as in Ref. ^. 
Figure ^ shows the homogeneous cavitation pressure as a function of T for the ^ values we 
have been using. 

In conclusion, we have developed a phenomenological model to size the effect of dis- 
sipation in the cavitation process in liquid ^He that allows one to handle realistic critical 
cavitation configurations near the spinodal line, and to treat both helium isotopes within 
the same frame, using the dissipationless limit of the method in the case of ^He. The results 
we have obtained indicate that for liquid ^He even a moderate dissipation may reduce the 
crossover temperature in a non- negligible amount, displacing the homogeneous cavitation 
pressure towards the spinodal value. Viscosity may then be the reason of the inconclusive 
results for quantum cavitation reported in Ref. |^ which, if confirmed, would indicate that 
dissipation plays a crucial role in quantum cavitation in liquid helium. The experimental 
study of cavitation in undersaturated ^He-^He mixtures might then uncover a structure much 
richer than that theoretically described in Ref. p, since ^He is still superfiuid and ^He is in 
the normal phase. This would open the possibility of studying the influence of dissipation 
in the cavitation process varying the ^He concentration. 
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TABLES 



TABLE L Equation of state, sound velocity, energy barrier and quantum action (^ 
the spinodal point. 



(A- 



p 

(bar) 



Cs 

(m/s) 



AJ7 

(K) 



0) near 



s/n 



.0123 
.0124 
.0125 
.0126 
.0127 
.0128 
.0129 
.0130 
.0131 
.0132 
.0133 
.0134 
.0135 
.0136 
.0137 
.0138 
.0139 
.0140 



-3.08 
-3.07 
-3.06 
-3.05 
-3.03 
-3.01 
-2.99 
-2.96 
-2.93 
-2.90 
-2.86 
-2.83 
-2.78 
-2.74 
-2.69 
-2.64 
-2.59 
-2.53 



42.3 

47.8 

52.7 

57.4 

61.8 

66.1 

70.1 

74.1 

77.9 

81.7 

85.3 

88.9 

92.4 

95.9 

99.3 

102.7 

106.1 

109.4 



1.3 

1.8 

2.5 

3.2 

4.0 

4.9 

5.8 

6.9 

8.0 

9.3 

10.6 

12.1 

13.7 

15.5 

17.4 

19.6 

21.9 

24.6 



16.0 

19.4 

22.9 

26.8 

31.2 

36.0 

41.5 

47.8 

55.0 

63.5 

73.5 

85.5 

100.3 

119.2 

144.0 

175.4 

212.6 

256.6 



FIGURES 

FIG. 1. T* as a function of pressure for different ^ values (in /iP). The homogeneous cavitation 
pressure Ph{T*) is shown as circles for {Vt)exp = 10® A s. 

FIG. 2. Effective quantum action in % units as a function of pressure for the values of ^ 
employed in Fig. 1. 

FIG. 3. Particle density profile po{r) in A " for P = —3 bar (solid line), as well as p^(r) 
densities for three ^ values (arbitrary units). 

FIG. 4. Homogeneous cavitation pressure Ph as a function of T for {Vt)exp = 10® A s, and 
the ^ values employed in Figs. 1 and 2. 
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